The breakdown of the Nagaoka phase in the 2D t — J model 
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In the limit of weak exchange, J, at low hole concentration, 5, the ground state of the 2D t — J 
model is believed to be ferromagnetic. We study the leading instability of this Nagaoka state, 
which emerges with increasing J. Both exact diagonalization of small clusters, and a semiclassical 
analytical calculation of larger systems show that above a certain critical value of the exchange, 
Jcr ~ tS^ , Nagaoka's state is unstable to phase separation. In a finite-size system a bubble of 
antiferromagnetic Mott insulator appears in the ground state above this threshold. The size of this 
bubble depends on 5 and scales as a power of the system size. A''. 



Recently, much interest was focused on the behavior 
of strongly correlated electron systems, which can not 
be explained by weak coupling perturbation theory. A 
variety of unusual phenomena such as, e.g. high tem- 
perature superconductivity and quantum magnetism, is 
believed to require a non perturbative descriptiontl. An 
important paradigm for the study of strongly interacting 
electrons in general, is the Hubbard model for interacting 
particles on a lattice. This is probably the simplest possi- 
ble model that captures some of the behavior of strongly 
correlated electrons. It is therefore widely used to study 
various correlation driven effects, which are not described 
by a perturbative approach. These include the metal- 
insulator (Mott) transition, and the superconductivity of 
the high-Tc compounds. 

The Hubbard model was originally introduced in an 
attempt to describe quantum ferromagnetism of itiner- 
ant electrons in narrow band metalsEl. However, it is now 
well known that it is also a model for quantum antiferro- 
magnetism. The effective Hamiltonian which governs the 
low-energy behavior of the Hubbard model for a nearly 
half-filled band in the U ^ oo limit is the t — J model of 
itinerant fermions on the lattice 
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The occupation number of each site rii — a\.^ai'\ H- al^^aii 
can be either (a hole) or 1 (a spin), since double- 
occupancy is forbidden by strong on-site Hubbard repul- 
sion. The spin operators Si are given in terms of the 
Pauli matrices aap 
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Since the spin exchange coupling constant J w At^ /U is 
positive, the on-site Hubbard interaction translates to an 
anti-ferromagnetic (AFM) super-exchange, which favors 
an AFM correlated ground state. 

Nevertheless, in one extreme case this model has a fer- 
romagnetic (FM) ground state, known as the Nagaoka 



stataj. A fully polarized state minimizes the total en- 
ergy of a single hole in an otherwise half-filled band, at 
least in the limit J ^ 0. This is due to the hole ki- 
netic energy which favors FM ordering. The nature of 
the ground state for finite J is therefore determined by 
the competition of this FM tendency with the AFM ex- 
change. 

In this work, we study the ground state of the 2D t — J 
model for low concentration of holes, as a function of 
the AFM interaction strength J . We aim to identify the 
leading instability of the Nagaoka state, i.e., the state 
which minimizes the energy for non- fully polarized states. 
Based on numerical evidence, we claim that the leading 
instability is towards the creation of an AFM bubble, 
while the holes are confined to the FM region. The tran- 
sition between Nagaoka-state and this "bubble" state is 
a first order transitionCj, including a jump AS" 3> 1 in 
the total spin. We will concentrate on the t — J model, 
Eq. (1). However, our study is focused on the small J 
regime, and thus our considerations apply equally to the 
Hubbard model. 

We start by summarizing in some detail what is known 
and conjectured about the nature of the t ~ J model 
ground state. For half-filling, where each site is occupied 
by one electron, no hopping is possible, and the Hamil- 
tonian reduces to the quantum Heisenberg Hamiltonian 
with AFM spin exchange. Although the 2D Heisenberg 
model can not be solved analytically, there is strong 
numerical evidence, (obtained from exact diagonaliza- 
tion and Quantum Monte Carlo simulations) that the 
ground statfe-kas long range AFM correlations at zero 
temperatureQEl. 

The presence of holes makes the picture more compli- 
cated. Each hole hopping creates changes in the spin con- 
figuration, unless the spin polarization is uniform. The 
resulting excitations in the system inhibit the hopping, 
since the hopping probability is reduced by a factor pro- 
portional to the overlap between the original and final 
spin wave functions. This effect results in narrowing the 
kinetic energy band, thus increasing the kinetic contri- 
bution to the energy. The band width is maximized in 
the fully polarized state for which the spin configuration 
is unaffected by permutations of different spins. Thus, 
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while the J term in the Hamiltonian favors an AFM or- 
dering, the kinetic t term favors a ferromagnetic (FM) 
state. The competition between these terms depends on 
the hole density and interaction strength J/t, and deter- 
mines the physical properties of the t — J model. 

The question of charge carrier dynamics on an AFM 
background was extensively studied using the string pic- 
ture, self-consistent Born approximation, and numerical 
studiesO. For J ^ t a single hole is unable to alter 
its AFM surrounding, and its kinetic energy comes only 
through coupling to the quantum fluctuations of the spin 
system. For moderately large J {J/t w 5 — 10) the hole 
develops a spin-polaron around itself. Within this po- 
laron, the AFM order is suppressed by hole hopping. The 
size of the polaron increases with decreasing J/t. The 
"string" picture emerges when this size is much larger 
than the lattice constant. In this picture, the hole acts 
as a particle with mass 1/t subject to an effective po- 
tential resulting from the string of flipped spins created 
by its hopping. This effective potential is a linear poten- 
tial with slope J. Accordingly, the size of the polaron 
scales like (J/t)~^/^. In this regime the hole movement 
scrambles its surrounding spins, thus suppressing the lo- 
cal AFM order parameter, but does not yet create a lo- 
cal niAffnetic moment larger than the spin- 1/2 of a single 
hole.BQ 

As soon as J/t decreases beyond some critical f^lue 
(Jc < 0.1), the nature of the ground state changegj. It 
turns out that lower energy is achieved by creating a 
FM polarona, where not only the AFM ordering is de- 
stroyed, but also FM correlations and a larger magnetic 
moment are formed in the vicinity of the hole. The ra- 
dius of the polaron in this regime scales like {J/t)~^/^. 
In the limit J ^ 0, the kinetic energy which favors FM 
ordering dominates the tendency towards AFM coupling 
of neighboring spins represented by the J term, and the 
size of this polaron diverges. As a result, the ground 
state becomes fully polarized. This statement was rigor- 
ously proven in the the celebrated paper of NagaokaO for 
a single hole case. It was shown that Nagaoka FM phase 
survives for sufficientUi,small density of holes, even in the 
thermodynamic limitll2l. 

Thus, this apparently simple model exhibits a 
ferromagnet-antiferromagnet transition at zero temper- 
ature, as the interaction strength J increases. Very little 
is known about the properties of the system in the transi- 
tion region. Even the dependences of the magnetization 
curve on J and the hole concentration S are yet to be 
determined. In the following, we try to understand some 
features of the transition. We focus on the boundary 
regime of Nagaoka's phase, i.e., the transition from full 
to partial polarization. 

Naturally, the question of the stability of the Nagaoka 
fully polarized state with respect to switching on the 
AFM exchange term, as well as finite hole density S, 
attracted much interest. Several variational wavefunc- 
tions were suggested in order to yield batrnds for the Na- 
gaoka stability region in S-U plane (Seetil and references 



therein) . Most of these estimates were based on the belief 
that tiie transition from Nagaoka state is continuous at 
T = dB. In other words, it was implicitly assumed that 
the leading instability of Nagaoka's state is a single spin 
flip (SSF). According to this picture, the transition to 
the AFM singlet ground state with increasing J/t occurs 
gradually, through small incremental AS* — 1 changes in 
the total spin. However, it was suggested that for suffi- 
ciently low S a phase separatiflp-instability, rather than 
a SSF one, might be relevantE3Ej. 

Many authors have discussed the possibility 
of phase .sepairaJiian in the Hubbard and t — J 
modelsl^JO'E^llSliZllla. For high concentration of 
holes (low electron density,), phase separation does occur 
for J sufficiently largalj'tj. However, the small J/t 
case is still under debate. Some groups argue that the 
ground state phase separates for all-jsalues of J/t for 
sufficiently low concentration of holealjO. Others claim 
that there exists a critical value Jps (estimates for its 
value vary between Jps OAt and Jps 1.4i), such 
that for J < Jps thCi-ground state is uniform even for 
vanishing hole densityli3. 

In the following, we present evidence that the leading 
instability for low densities is indeed a phase separated 
state. We therefore claim that phase separation does oc- 
cur for small values of J/t, given that 6 is sufficiently 
small. We further show that for a finite-size sample the 
transition in the ground state between the Nagajaka phase 
and the phase separated state is discontinuousllj, includ- 
ing a large change in the total spin AS* >> 1. Thus this 
transition can not be captured by the SSF variational 
studies. It is therefore reasonable to assume that better 
bounds for Nagaoka stability might require considering 
many simultaneous, rather than sequential, spin-flips. 

We start by examining at the single hole case. It is 
straight forward to show that Jcr, the J value needed to 
destabilize Nagaoka's jState, for two spin flips is smaller 
than for one spin flipllj. This can be demonstrated by 
deriving an effective Hamiltonian of a single hole and 
flipped spins, in the background on the FM Nagaoka 
state. Since the number of flipped spins is presumably 
small, it is convenient to describe the system in terms of 
the two types of "particles" - the hole and the flipped 
spins. The kinetic part of the Hamiltonian induces hop- 
ping of the hole to its nearest neighbors, with an effective 
mass l/{2t). The non diagonal part of the exchange in- 
teraction induces hopping of the flipped spins to their 
nearest neighbors, with an effective mass 1/J. The ratio 
between the effective masses of the spin and the hole is 
therefore J/2t, which is small around Jcr- Therefore, one 
can assume that the spin is static for calculating its effect 
on the hole energy, i.e., use Born-Oppenheimer-like adia- 
batic approximation. The diagonal part of the exchange 
term contributes the energy — J/2 per flipped spin. In 
addition to these terms, one should take into account the 
constraint that each site can be occupied by no more than 
one "particle" (hole or flipped spin). For that purpose, 
one can introduce effective, inflnitely strong, on-site re- 
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pulsion between the particles. There are other terms in 
the effective Hamiltonian resulting from nearest neigh- 
bors interaction between the particles, but they do not 
contribute to the energy in the leading order in J /t. The 
Nagaoka state is destabilized as soon as the 0( J) mag- 
netic energy gained by the flipped spins overcomes this 
kinetic energy gain. 

For a state with two spin flips, the kinetic energy in- 
crease depends on the distance between the spin flips. 
When they are far apart, the kinetic energy lost by the 
excluded area is about twice the value for a single flip. 
However, as the two spin flips get closer, the effective 
excluded area for the hole decreases. When the distance 
is small compared to the hole wavelength, the excluded 
area is not changed much, and the kinetic energy increase 
is about the same as for one flip. Thus, there is an effec- 
tive attraction between the flipped spins, and they bind 
together, with total energy similar to that of a single 
flip. On the other hand, the magnetic energy is simply 
proportional to the number of spin flips, and does not 
depend on the distance between the spin flips, as long 
as they are not nearest neighbors. Thus the magnetic 
energy for two flipped spins, bound or not, is twice as 
much as the energy for a single flip. Therefore, Jcr, for 
which the magnetic energy balances the kinetic energy 
increase, is lower for a state with two flipped spins. One 
concludes that, for a single hole, the transition between 
the maximum spin Nagaoka state and a lower spin state 
involves a spin jump A5 > 1. 

In order to describe the leading instability of the Na- 
gaoka ferromagnet with respect to switching on the J 
term one should, therefore, determine the number of spin 
flips which minimizes Jcr. We start with presenting re- 
sults of exact diagonalization for small rectangular axb 
clusters. We choose periodic boundary conditions when 
the large axis of the cluster has even number of sites 
(a — 2n), and anti-periodic boundary conditions other- 
wise (a = 2n -I- 1). Under these conditions the J = 
ground state is always fully polarized. For each of the 
following clusters we diagonalize exactly the full many- 
particle Hamiltonian, for different spin sectors and var- 
ious J values. We determine Jcr and the value of the 
grojund state spin for J = Jcr+- It was already pointed 
outo, based on exact diagonalization results, that for big 
enough clusters, large number of spins may flip together. 
We present here a more systematic study of finite clus- 
ters, showing a large spin jump AS* > 1. For the one hole 
case, the biggest cluster studied was a 6 x 4 torus. The 
effective size of the full Hilbert space can be substan- 
tially reduced by excluding the doubly occupied states 
and using spin and translational symmetry. For 6x4 
torus, this reduces the problem to the diagonalization of 
a 1352078 x 1352078 matrix. We employed the Lanczos 
algorithm to carry out this diagonalization. The results 
for the different clusters are summarized in table I. Ex- 
tension of these exact diagonalization studies to larger 
systems or more than one hole is limited by the size of 
the Hilbert space. 



Cluster 


AS 


J CT 


Cluster 


AS 


J CT 


2x8 


3 


0.0225 


4x4 


4 


0.0629 


2 X 10 


4 


0.0134 


4x5 


6 


0.0398 


2 X 12 


5 


0.0086 


4x6 


6 


0.0269 



TABLE I: Exact diagonalization results for the breakdown on 
Nagaoka's state for the case of one hole. AS is the jump in 
the ground state spin at the transition. 



We did repeat the calculation for two holes on the 
above clusters. It is well known that for two holes,|Oii 
a torus the ground state is a singlet even for J = 0tBE3. 
This is due to the fact that a slowly varying, locally polar- 
ized, spin background creates the effect of a fictitioiisJlux, 
which minimizes the kinetic energy of the two hole£j. As 
it was done for one hole, we chose the boundary condi- 
tions to overcome this effect. However, in all cases stud- 
ied, a ground state with small nonzero spin was never 
observed. We always dealt with an abrupt transition to 
the singlet S — state, AS* being as large as the largest 
possible spin, S'max, even larger than the jump in the 
spin for the one hole case. This, however, is a particular 
feature of systems with two holes, which might not be 
relevant for the behavior of large systems. 

These exact results show that the number of spins 
which flip at the transition is significantly larger than 
one, and it increases with the system size. They indicate 
that the breakdown of Nagaoka phase might indeed oc- 
cur through an abrupt spin change. It looks interesting 
to understand how A^ scales with system size. Besides, 
it is not clear to what extent are the results specific for a 
particular small number of holes (1 or 2), i.e., extension 
of the results to more holes is required. Clearly, exact 
diagonalization is not an adequate tool to clarify these 
questions. 

To get an insight into the behavior of larger systems, 
we use the spin-hole coherent state path integi'al formal- 
ism. We follow the derivation presented inE3. A semi- 
classical approximation for the ground state is obtained 
through the formal large S expansion. In the classical 
limit, S oo, the spins are frozen (i.e., the spin part 
of the coherent state path integral is time-independent). 
The energy is determined solely by the interplay between 
hole dynamics and the classical AFM interaction between 
the spins. In essence, this limit deals with a classical 
spin field, interacting with the quantum holes. Clearly 
this approximation can not capture the full complexity of 
the exact wave function, as it ignores dynamical correc- 
tions to the spin background as introduced by, e.g., the 
dressed holes. Thus, some of the exotic phases suggested 
for the t — J model are essentially not taken into account 
in this approximation. However, in some cases the insta- 
bility towards these phases can be described within the 
semi classical approximation. For .example, the incom- 
mensurate AF, or spiral phaseHEJ, with dressed holes, 
supercedes the Nagaoka phase for two holes. Within the 
semi classical approximation, ignoring the dressing of the 
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FIG. 1: The spin configuration of tfie leading instability as 
obtained from a semi-classical calculation. The hole density 
is 25/256. 16 spins flip in a slanted square configuration. The 
hole density vanishes inside the AFM bubble, while being 
approximately constant in the FM region. 

holes, the incomensurate AF background becomes degen- 
erate with the Nagaoka phase. Thus one can expect to 
find signs for the emergence of complex spin structures 
even within this simplified approach. The same approach 
was previously used to study the formation of the FM 
polaron. Here we use the same technique to study the 
extreme low J case J = Jci+, at the breakdown of Na- 
gaoka phase, where the FM region occupies most of the 
system. We use a Monte-Carlo algorithm to find the spin 
configuration which minimizes the sum of the energy of 
the holes and the magnetic energy of the (classical) spin 
field. 

We study different lattice sizes (up to 16 x 16) and 
various number of electrons close to half-filling {0 < 6 < 
0.12). In all cases, we find the same behavior: (a) all 
the spins align co-linearly (b) we do not see any signa- 
ture suggesting the emergence of exotic spin configura- 
tions (spiral, canted, etc.) for these small values of J, 
and most importantly (c) the uniform Nagaoka's state 
breaks down into a phase separated state, with a hole- 
rich FM region and an AFM bubble with no holes inside. 
The size of the bubble at the transition is large relative 
to the lattice constant, giving a big jump in the spin. As 
an example we present in figure |l| results for a 16 x 16 
lattice with 25 holes. These results, for different lattice 
sizes and densities, suggest that Nagaoka state breaks 
down by forming an AFM bubble, whose size is large 
compared to the lattice constant. Contradicting evidence 
were recently presented ijpgarding the existence of phase 
separation for small JtZlEEl. We find that the above de- 
scribed phase separated state seems to be stable even for 
J << t, provided that the haLcj-dcnsity S is sufficiently 
small, in agreement with Ref.EjO. The stability of the 
phase separated state for larger J, and its relation to the 
above discussed phase, is beyond the scope of this study. 

Motivated by these results, we turn to study the sta- 
bility of the Nagaoka phase with respect to phase sepa- 
ration. That is, we compare the energy of the Nagaoka 
state with that of a system with an AFM ordered do- 



4 

main. To first approximation, the AFM domain acts as 
an infinite potential barrier, i.e. the holes are confined 
to the FM region. Due to this confinement, the kinetic 
energy of the holes increases with the size of the AFM do- 
main. On the other hand, the AFM domain contributes 
a magnetic term to the energy, which also increases with 
its size. We estimate the magnetic energy per site in the 
AFM domain by the value obtained in an infinite Heisen- 
berg system, 

E^,^ = -aJ{S + V/2), (3) 

where S is the area (number of sites) of the AFM do- 
main, V is its perimeter, and a is the energy per site of 
the ground state of the 2D Heisenberg model (Spin wave 
theory gives a ~ 1.1705, jvhich is in a good agreement 
with the numerical resultn) a = 1.169. An estimate to 
the increase in the holes' kinetic energy for a large system 
is obtained through the well known Weyl formula for the 
density of states. The number of levels up to an energy 
E for a free particle in a continuous 2D domain with area 
A and boundary perimeter V is given by 

{N{E))^1-{AE-V^ + IC) , (4) 

where Dirichlet boundary conditions are assumed, and K. 
is a constant term containing information on the geome- 
try and topology of the domaiio. The units are chosen 
such that h = 2m = 1 . Equivalently, for the Hubbard (or 
t~ J) model, the same formula holds near the bottom of 
the band, where the energy is in units of t and the bottom 
of the band is taken as the energy zero. In the following 
t is taken as the unit of energy, and the lattice constant 
is taken as the unit of length. The average {p{E)) of the 
density of states (DOS) p{E) = dN{E) jdE, is thus given 

by 

(P(^)> - ^ - (5) 

In the Nagaoka FM state, the available domain for the 
holes is the whole torus A = N . The DOS is therefore 
just the familiar Aj ^'n term and the energy is given by 

i^Nag = 2'kN6'^ = 'i'^Nl/N , (6) 

where Nh is the number of holes and TV is the total num- 
ber of sites. As soon as the AFM bubble is formed, the 
energy increases due to two reasons. First of all, the 
available area is reduced to A — N — and therefore 
the denominator in Eq. (^ decreases. Due to this fact, 
the total energy increases by the factor N / [N — S) . An- 
other contribution comes from the boundary term which 
reduces the DOS even further. Bearing in mind to com- 
pare the energy increase with the magnetic energy which 
is 0{S), we realize that the increase due to the boundary 
term 0{V) divided by S becomes singular (~ S^^/"^) for 
small S. We will see below that the transition occurs for 
S values much smaller than the system size. As a result. 
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FIG. 2: An illustration of the total energy (^) for different 
values of j. (a) No meta-stable state, (b) A meta-stable state 
appears . (c) The transition point, (d) J > Ja. 



boundary terms are important even for large N. The 
topological K. term also changes due to the AFM barrier. 
However, it can be checked explicitly that it does not 
affect the — > oo asymptotic behavior. 

Using the above averaged DOS, one can calculate 
the total energy of holes as a function of their density. 
Adding the magnetic energy, the total energy relative to 
the Nagaoka energy is given by 



E{S,s) = 2TrS^N 



2p 



-aJN s + 



1 - s 3yMN{l - s)3/2 
P 



2VN 



where we introduced the normalized quantities 



p = V/VN, s = S/N. 



(7) 



(8) 



It follows from the dimensional analysis that p = Cy^, 
where C is a dimensionless constant depending on the 
shape of the bubble. Although we deal with a lattice 
problem, in the limit we are considering of low doping and 
large bubbles, the perimeter term from the hole kinetic 
energy dominates. Since the dispersion relation for low- 
energy holes is asymptotically isotropic, this produces an 
isotropic perimeter energy so the lowest energy bubbles 
are circular (C — Zy^). It is straight forward to see that 
given the grand canonical potential, G = E{S, s) — N^S, 
with a fixed chemical potential fi, that the system goes 
through a first order phase transition. This transition 
occurs at Ja = /i^/(87ra), where the density jumps from 
S- = ///(47r) at J = JcT— to (5+ = at J = Jcr+- In 
the same time s jumps from s — (Nagaoka phase) be- 
low the transition, to s = 1 (antiferromagnet) above the 
transition. If the density is instead held between these 
two values, the ground state, ignoring surface terms, is 



a phase separated state, which is a mixture of the two 
phases. In the following we study the formation of the 
AFM bubble in a finite system, given the surface terms. 
The energy (0) as a function of s is presented in figure ^ 
for various values of J. At the origin, E{0) = 0, and the 
energy increases with s. For sufficiently large J, the func- 
tion has a maximum at which the energy is positive and 
a (local) minimum at s = s* (see Fig. ^ curve (b)). When 
J increases beyond a critical value J > Ja-, the value of 
the energy at this minimum becomes negative, the lo- 
cal minimum at becomes the global minimum of the 
function. The function crosses the x-axis at s — sq < s^, 
(see Fig. || curve (d)). For s > s* the function again 
increases. Therefore, as long as J < Jcr, the global en- 
ergy minimum is obtained at s = 0, corresponding to 
Nagaoka FM state, although a meta-stable state with an 
AFM bubble does exist for a range of J near, but be- 
low Jcr. As J increases beyond Jcr, the global minimum 
shifts to s*. It is convenient to introduce the variable j 



] 



J = 27r5Va(l+j)^ Jo"(l+.?-) (9) 



Jq' is the value of Jcr for an infinite size system. Thus, 
j measures the relative distance of J from the value of 
Jcr obtained ignoring boundary terms. In terms of this 
variable, one can express the large N asymptotics of sq 
(the value of s where the energy vanishes) and s», which 
minimizes the energy. 



So 



CiA/j -ly/m; A = 4/(3V7r(5) - 1 (10) 



1 



1 



(11) 



note that A depends only on the hole density 6. Thus, 
for J > 0, the minimal size of the droplet needed to 
destabilize Nagaoka's state remains finite in the N — > 
oo limit (5'o = SqN = 0(1)), while the global min- 
imum is obtained when the droplet is of macroscopic 
size (5, = s^N = 0{N)). Comparing these two expres- 
sions, one immediately sees that the transition occurs for 
j - 0(iV-i/3). 

In this regime sq and s* are connected with j via the 
equations 



j = 2s, + 



CA 



So 



CA 



(12) 



(13) 



The last equation has a solution provided that j > jcr = 
3{C'^A'^/16Ny/^. The transition point Jcr is therefore 



Jc: 



271,52 



l + 3((C^)Vl6A^)^/^ 



(14) 



Beyond the transition point J > Jcr, and for some s the 
energy of the bubble is negative and therefore the FM 
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FIG. 3: The area of the AFM bubble at the transition point, 
Scr, normalized by system size, as a function of A'^ for various 
N values. The solid lines are the asymptotic expansion up to 
order N~'^ for the same five densities (top to bottom) . 
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FIG. 4: Schematic phase diagram in the vicinity of the break- 
down of Nagaoka phase. Region I - Nagaoka phase is sta- 
ble. Region II - Bubbles bigger than the critical value are 
metastable. Region III - Bubbles bigger than the critical value 
are stable. 



state is unstable. At the transition point, the area of the 
optimal AFM bubble is 

So = Nso = Ns^ = iCAN/4f/^ (15) 

For j < jcr Eq. ( p^ ) has no solution and the energy is 
minimized by Nagaoka state. However, the energy func- 
tion does have a local minimum at for j > j^s = 
2~^^^jcr, and thus a large meta-stable bubble can be cre- 
ated. The size of the meta-stable bubble at j — jms+ is 

= s,N = (CAiV/16)2/3 (16) 

In fact, it can be shown that the N'^/'^ power law fol- 
lows from a very general argument. Whenever phase 
separation occurs between two phases, with densities 6i 
and ^2, the fraction of the first phase in the phase sepa- 
rated ground state of an infinite system at a fixed density 
Si < S < S2 is given by (^2 ^ S)/{S2 — Si). Thus, for den- 
sities very close to S2 the size of the bubble of the first 
phase is arbitrarily small. However, the surface energy 
term raises the energy of the phase separated state as 
compared to the uniform state, and thus finite size cor- 
rections arise. The energy cost to create a boundary be- 
tween the two phases, is proportional to (in d dimensions) 
yid-i)/d^ where V is the size of the bubble. This en- 
ergy makes it favorable to retain the uniform phase even 
at densities slightly smaller than S2- The formation of 
a bubble becomes energetically favorable when the den- 
sity is shifted away from the (infinite system) transition 
point, such that the difference in the volume energy of the 
uniform state as compared to the phase-separated state 
overcomes the surface term. In the generic case, this dif- 
ference is proportional to V'^/N where N is the system 
size, and V is the size of the optimal bubble. Compar- 
ing the two energies, one observes that at the transition 
V ~ Ard/(d+i). 



We performed a numerical study to validate the above 
calculations. The single-particle spectrum of a tight 
binding square lattice model with periodic boundary con- 
ditions and an excluded domain, was calculated for dif- 
ferent sizes of the excluded domain. We then calculated 
the ground state kinetic energy (which is the sum of the 
lowest TO eigenvalues, where m is the number of holes), 
and compared the increase of kinetic energy as a result 
of the excluded bubble with the gain in magnetic energy. 
For each size of excluded domain, a minimal value of J 
is needed to balance the increase of kinetic energy, thus 
stabilizing the AFM bubble. Jcr is determined as the 
minimum of these J values, i.e., the lowest value of J 
which allows for a stable bubble. The size of the bubble 
at this minimum is Sq. The area of the stable bubble at 
J = Jcr+ is presented in figure ^ as a function of S and 
N. It turns out that the leading asymptotic dominates 
only for a very big system, beyond numerical capabili- 
ties. We therefore calculated the sub-leading corrections 
0{N~'^/'^) and 0{N~^) and compared the numerical re- 
sults with this asymptotic expansion. The agreement 
with the analytical result is quite good for sufficiently 
large lattices, and it seems that the results converge to 
the asymptotic estimate as the system size increases, in- 
dicating that corrections due to deviations from the free 
particle band shape are not important for S < 0.05. 

The picture emerging from these calculation is as fol- 
lows (see also Fig |). At J = Jcr (Eq. ([l4|)) the Nagaoka 
state breaks down and one large AFM bubble with the 
area of order of 0(iV^/^) is formed (Eq. (|l5|)). As the 
size of the bubble scales only sub-linearily with system 
size, the magnetization per site is continuous in the ther- 
modynamic limit. In the vicinity of the transition point, 
only these large bubbles are stable, and therefore a large 
fluctuation is needed to destroy the Nagaoka phase. It is 
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therefore a meta-stable phase. For larger J, sq, the crit- 
ical size for a stable bubble, decreases, and for j ~ 0(1) 
(i.e. J — Jcr ~ Jcr) it becomes finite and system-size 
independent. This behavior is typical for a first order 
phase transition, where a domain of a new phase has to 
be large enough to survive. However, the jump in the 
magnetization is not extensive. 

Why doesn't the system create many more bubbles, 
once the first bubble is formed, thus reducing the mag- 
netization even further ? The reason is, as can be shown 
in Eq. (^, that the kinetic energy is not a linear function 
of s. Therefore, the energy cost to create two bubbles is 
more than twice the cost to create one. On the other 
hand, the magnetic energy gain, is linear in s, and there- 
fore it would just double by the creation of the second 
bubble. 

There is an ongoing discussion in the literature regard- 
ing the existence of a striped ground state in the Hubbard 
and t ~ J models. Experimentally, there is evidence for 
stripe modulations in doped cuprate^3, which are gen- 
erally believed to be described by the t — J model. Some 
authors found that the t — J model grownd state is in- 
deed striped for a wide range of dopingSa, while others 
claim that uniform or phase separated states have lower 
energyBa. According to the latter view, the origin of the 
experimental observation is attributed to the competi- 
tion between the local tendency for phase separation and 
the long-range Coulomb interaction, which is missing in 
the t— J modelEZl. 

In this study, we compared the energy of a striped state 
with the the energy of a phase-separated state with an 



AFM bubble. Within our approach, the kinetic energy 
increases in the striped phase due to the higher surface 
energy, with nothing else to compensate for this increase. 
Accordingly, the striped state energy is higher, suggest- 
ing that the long-range Coulomb repulsion is needed in 
order to create a striped ground state for J/t <C 1. Our 
approach can not rule out the possibility of a striped 
ground state of the t — J model aX J ^ t. 

In conclusion, we presented analytical arguments, ex- 
act diagonalization results, and semi-classical calcula- 
tions of the 2D J model, which suggest that at small 
hole concentration, 5, and rather weak AFM coupling, 
J, the Nagaoka ferromagnetic state becomes unstable to- 
wards the creation of an AFM bubble. In this phase sep- 
arated state, the holes are confined to the FM regime. 
At the transition only a single large 0{N'^^^) bubble is 
(meta)stable. Thus, the magnetization is continuous at 
the transition in the thermodynamic limit. However, the 
jump in magnetization per unit area in a system with fi- 
nite number of sites iV, scales as N'^^'^. This dependence 
of the critical bubble size, and thus of the magnetization, 
on the size of the system is a typical finite size effect in 
a phase separated ground state. 
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